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Quenched lattice QCD calculations with Wilson-type fermions at small quark masses are impeded by exceptional 
configurations. We show how this problem can be resolved in a practicable way without changing the physics in 
the chiral limit. 



1. INTRODUCTION 

We consider Wilson fermions with or without 
0(a) improvement. The fermionic action reads 

S F = K = W + M ' f 1 ) 

where W is the Wilson-Dirac operator: 

W =fl + X (2) 

with 

1 



D» = ^{D+ + D-l 
X = --D+D- + ^c sw <V,iV, 



(3) 
(4) 



£)+ (D^ ) being the forward (backward) lattice 
derivative, and where AI is the Wilson mass term: 



M = — - 4. (5) 
2k 

The critical value of the hopping parameter, k c , 
is taken to be the limiting value at which the en- 
semble averaged pion mass vanishes. We call this 
limit the chiral limit. 
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*Poster presented by G. Schierholz at Lattice 1998. 



10 20 30 40 50 60 
t 

Figure 1. The pion correlator on the 32 3 64 lat- 
tice at fi = 6.2 and k = 0.1354 for improved 
fermions with csw = 1.614. At these parame- 
ters m^/nip ~ 0.5 which corresponds to a quark 
mass of w 40MeV. The solid symbols (and line) 
represent the ensemble average with the excep- 
tional configurations omitted. The dashed line 
is the result of a single, particularly exceptional 
configuration. 
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Figure 2. A sketch of the eigenvalue spectrum of W c (before the rotation) and W<f w (after the rotation). 



In quenched lattice calculations with Wilson 
fcrmions one encounters exceptional configura- 
tions [jjj which give rise to poles in fermionic ob- 
servables, such as hadron correlators, at small 
quark masses. The origin of these poles are real 
eigenvalues Aj of the Wilson-Dirac operator W at 

A,<4-^-. (6) 

Exceptional configurations are particularly aggra- 
vating for improved fermions, and they have lim- 
ited calculations to larger quark masses and larger 
values of (3. We expect to find such poles also in 
the full theory, i.e. with dynamical fermions, at 
hopping parameters n ^ n S ca, where K soa refers 
to the sea quark mass. 

In Fig. 1 we give an example of an exceptional 
configuration which we encountered in a recent 
simulation with improved fermions We see a 
dramatic deflection from the anticipated ensem- 
ble average. Similar peaks are observed in the 
p and nucleon correlators. It is quite clear from 
this figure that the problem of exceptional config- 
urations cannot be solved by accumulating larger 
statistical samples. 

Bardeen et al. (|] identify the eigenvalues @ 
with the zero modes of the continuum theory, 
and they propose to shift each such Ai to = 
4 — 1/2k c . There are doubts, however, that this 
identification is correct in general ||]. An ele- 
gant solution would be to employ Neuberger's ac- 
tion H . But whether this is numerically feasible 



has still to be seen. 



2. PROPOSAL 



The solution we propose makes use of a free- 
dom in regularizing lattice fermions. It has the 
advantage that it is easy to implement numeri- 
cally, and that it does not change the physics in 
the chiral limit so that the improvement program 
a la Sheikholeslami and Wohlert |(| remains un- 
altered. 

We rewrite K. in (Q) as 



m, 



(7) 



where W c and m are the critical Wilson-Dirac 
operator and the physical bare mass, respectively: 

1 

2k, c 

i i 



w c =y> + x c , x c = x 



-4, 



(8) 
(9) 



We then apply a chiral rotation to W c , 

W c -> W* = e^Wj^ , 

= $> + X c (cos 2^ + 175 sin 20), (10) 

and take JC* = W^+m to be the new fermion ma- 
trix. The resulting action is equally well founded 
as (|l|) and (0) Q], though, in general, it breaks 
parity invariance. In Fig. 2 we show the effect of 
the transformation on the eigenvalue spectrum of 
the Wilson-Dirac operator for <j> = 7r/4. In this 
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Figure 3. The pion correlator for a single ex- 
ceptional configuration on the 16 3 32 lattice at 
f3 = 5.7 and c sw = 2.25, k = 0.1295 before 
(o) and after the rotation by 4> — ir/4 (•), with 
k c = 0.13074 The data for improved fermions 
are compared with the result for Wilson fermions 
(O) at an equivalent quark mass. 



case W c is anti-hcrmitian, and all eigenvalues are 
mapped onto the imaginary axis, as is the case 
for staggered fermions. 

3. TEST 



The rotation (10) is mathematically equivalent 



to the transformation of the mass: 
m — > = m (cos 2<j) — 1^5 sin 20). 



(11) 



In the following we shall use this transformation. 
Furthermore we take 4> — 7r/4. But we found that 
the method works equally well for angles as small 
as \(f>\ = 0.05. 

To test the idea, we have looked at exceptional 
configurations which we encountered in our runs 
with improved fermions || on the 16 3 32 lattice at 
(3 — 5.7. If the method works at this coupling, we 
expect that it will work also at larger values of f3. 
In Fig. 3 we show the pion correlator for a single 
exceptional configuration before and after the ro- 
tation. Before the rotation we find a similar peak 
in the pion correlator as seen in Fig. 1. After the 



rotation we find a correlator which is perfectly 
well behaved. It turns out to be in good agree- 
ment with the ensemble average ||. Note that 
in first approximation the physical quark mass 
is left unchanged by the rotation. For compar- 
ison we have computed the pion correlator also 
for Wilson fermions. The hopping parameter was 
chosen so as to give roughly the same pion mass. 
We find good agreement with the rotated action. 

4. DISCUSSION 

From both, Figs. 1 and 3 it follows that the 
(suspected) near-zero eigenmode of K, acts as a 
source term, and that the excitation it induces 
decays exponentially with roughly the mass of the 
pion. We would not expect to find such a behavior 
for an instanton configuration. There is also no 
sign of an abnormal fluctuation found in the case 
of Wilson fermions. 

We conclude that the problem of exceptional 
configurations may be resolved by a chiral rota- 
tion of the critical Wilson-Dirac operator. The 
method allows calculations at any quark mass. 
On top of that, the results of the improvement 
program remain valid. 
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